1 



Progress of Theoretical Physics, Vol. **, No. **, January 2002 

Dynamical Effects beyond Mean Field in Drip Line Nuclei 

Hiroyuki Sagawa*) Gianluca Colo 1 and Toshio Suzuki 2 

Center for Mathematical Sciences, University of Aizu, 
Aizu-Wakamatsu, Fukushima 965-8580, Japan 
1 Dipartimento di Fisica, Universita degli Studi and INFN 
Sezione di Milano, Via Celoria 16, 20133 Milano, Italy 
2 Department of Physics, College of Humanities and Sciences, Nihon University 
£NJ \ Sakurajosui 3-25-40, Setagaya-ku, Tokyo 156, Japan 

' (Received January 11, 2002) 



(N 
(N 



New structure problems are raised in nuclei near drip lines by recent RIB experiments. 
We pointed out the importance of dynamical effect beyond mean field theories to the new 
structure problems. 



. 81. Introduction 

Recently, nuclear structure in neutron-rich nuclei has attracted much attention 
because of its exotic nature, such as halo and skin 1 )" 3 ). Recent advances of sec- 
| ondary beam techniques allow us to measure useful structure information on masses, 

interaction cross sections and excited states. The shell structure is one of essential 
issues in nuclear structure. Recently, the magic numbers in the neutron-rich region 
have been extensively studied both experimentally and theoretically. For example, 
the melting of the N=8 closed shell was pointed out due to the large mixing of lpi/2 
and 2si/2 orbitals in n Be 4 ), n Li 5 ) and 12 Be 6 ). The disappearance of the N=20 
closed shell was also shown in 32 Mg experimentally 7 ) . Very recently a possible new 
magic number is pointed out at N=16 by the analysis of the neutron separation 
^ ' energies and the interaction cross sections 8 - 1 . 

So far the nuclear structure problems in light nuclei have been studied mostly by 
using mean field theory 9 ) and shell model calculations 10 ) . It is known that the corre- 
lations beyond the mean field have substantial effects on the single-particle energies 
and electromagnetic transitions. In particular, the effect of collective excitations has 
been studied by the particle- vibration coupling and found to be important 11 » 12 ). In 
this paper, we present the study of dynamical correlations on single-particle energies 
in light neutron-rich nuclei based on the 10 Be and 24 core by using the particle- 
vibration coupling model 13 ) . The polarization charges for electric transitions in drip 
line nuclei have been discussed also by using the particle-vibration coupling model 
elsewhere 14 )' 15 ^. In the present version of the model, firstly Hartree-Fock (HF) cal- 
culations with Skyrme interactions are performed to obtain the zeroth-order single- 
particle energies and the wave functions. Secondly, we couple to them the natural 
parity vibrational states with J 7r =l _ , 2 + and 3 _ by using the particle-vibration 
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coupling Hamiltonian derived from the same Skyrme interaction, and the HF wave 
functions. In Section 2, we describe our model. Numerical results are given in Section 
3. Summary is given in Section 4. 

§2. Particle-vibration coupling model 



(a) 



(b) 



The particle- vibration coupling Hamil- 
tonian has the form 

R pv = E EH«rMr)F AM (f)|/3) at 0/J , 

a/3 Xnfi 

(2-1) 

where a, (3 label the single-particle 
states, v(r) is the form factor asso- 
ciated to the particle-hole interaction 
V p h derived from the Skyrme force 
(V p h{ri,r 2 ) = v(ri)6(ri -r 2 ) 16) ), and 
the radial transition density 5g n (r) can 
be derived from the RPA phonon. 

The second-order correction due to 
the RPA phonons to the single-particle 
energies coming from the perturbing 
Hamiltonian (2T) is obtained by the real 
part of the four diagrams depicted in 
Fig. 1. The imaginary part, which gives 
rise to the width of the single-particle 
states, is not discussed in the present 
work. The corresponding energy shifts 

for the diagrams can be found, e.g., in Ref. 17 \ and are given to be, 



(c) 



(d) 



Fig. 1. The four diagrams corresponding to 
the second-order couplings of particles (left 
column) or holes (right column) with vibra- 
tions, represented by wavy lines. See the 
text for details. 



Ae a = ReJ2 



\,n 



2ja + 



^ \Sdru a (r)u (rMr)S 0n (r)\ 2 \(a\\Y x \\l3)\ 2 | 
1 f3>F £ 



E 



£a ~ (£/3 + Wn) + ^ 

\fdru a (r) U p(r)v(r)5 0n (r)\ 2 \(a\\Y x m\ 2 . 



, (2-2) 



where e a and £3 are HF single-particle energies, u a and up are the radial HF wave 
functions and u> n is the RPA phonon energy. The sum over (3 spans unoccupied 
particle states (occupied hole states) when the notation j3 > F (j3 < F) is employed. 
If the state a is a particle state, the first and second terms in the above equation 
correspond to the upper and lower left-hand side diagrams (a) and (b) in Fig. 1. If a 
is a hole state, the first and second terms correspond, respectively, to the lower and 
upper right-hand side diagrams (d) and (c) in Fig. 1. On quite general grounds, it 
is expected that the real parts of Eq. (2-2) are larger for the "leading" diagrams (a) 
and (c) shown in Fig. 1 than the diagrams (b) and (d) in Fig. 1. This is because the 
energy denominators of the upper row diagrams (a) and (c) are expected to be smaller 
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than those of the diagrams (b) and (d). Furthermore, the energy denominator of the 
contributions (b) and (d) have opposite sign, so their contributions should cancel 
partially those of the leading terms. In 24 we will find exceptions to these rules due 
to the possibility that the leading contributions are not possible because of phase 
space limitations (hole intermediate states may not be available in these light nuclei 
if large angular momenta are required by selection rules). This will be one of the 
important qualitative messages of the present work. 



§3. Results 



We present numerical results for 
the single-particle energies in the 10 Be 
and 24 core. The shell structure 
of 10 Be core in the present study is 
compared with the results of Ref. 
(see also the coupled-channel calcula- 
tion of Ref. 19 ) ) . Since it is somewhat 
questionable to apply the self-consistent 
HF+RPA approach to obtain realistic 
single-particle energies and collective ex- 
cited states in very light nuclei, we adopt 
a phenomenological approach for 10 Bc. 
Namely, we start from a Woods-Saxon 
potential as the zeroth-order potential 
and calculate the particle-vibrational 
coupling with empirical phonon energies 
and transition strengths. 

In 10 Be, the existing Skyrme forces 
give a HF single-particle spectrum with 
a large IP1/2-IP3/2 g a P- The interaction 
SIII gives, for instance, e(lp 1 / 2 )=-3.80 
MeV and e(lp 3/2 )=-11.13 MeV. Using 
a standard Woods-Saxon parameteriza- 
tion 20 ), 

U = (-V +V T ^^)f(r)+VUl-s)R 2 ^- 

(3-1) 

(with Vq = 55 MeV, V T = 33 MeV, VJ S = 
ai s V and cti s = 0.44, and where f(r) is 
the usual Fermi function with ro = 1.27 
fm and a = 0.67 fm), the energy gap 
is reduced but the Fermi surface lies in 
a similar position as in the Skyrme-HF 
case. In fact, the Woods-Saxon results are 
MeV. This lpi/2 single-particle state, togel 
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Fig. 2. The neutron single-particle spectrum 
in 1 Be core calculated by a Woods-Saxon 
potential plus particle-vibration coupling. 
The energies are obtained, following the 
levels from left to right, by the Woods- 
Saxon potential, and including the effect of 
the 2]*~ state, of the and 2^ states, and 
of the 2+, 2^~ and 3J - states, respectively. 
The excitation energies and B(E2) values 
of the 2+ and 2^~ states are taken from ex- 
perimental data 21 ', while those of the 3^ 
state are obtained by shell model calcula- 
tions. The cutoff radius for the calcula- 
tion of the single-particle wave functions is 
taken to be 20 fm. Values corresponding 
to 15 fm are given in parenthesis for the 
sake of comparison. See the text for more 
details. 

e(lp 1/2 )=-3.93 MeV and e(lp 3/2 )=-9.36 
her with the 2si/ 2 unoccupied state, are 
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shown in the first column of Fig. 2. 

The vibrational states included in our model space are limited to those for 
which experimental information is available. The energy and reduced transition 
probability B(E2,0 — > 2^) of the lowest 2f state in 10 Be has been measured to 
be 3.37 MeV and 52 e 2 fm 4 , respectively 22 ) . The effect of coupling between the 
single-particle states and this phonon can be seen in the second column of Fig. 
2. The coupling reduces the 2s 1( / 2 -lpi/2 gap from 5 to 2 MeV, melting the N=8 
shell closure. Interpreting those as the ^ + and \ levels of 11 Be, the shifts in- 
duced by phonon coupling are not enough to give an inversion of the two levels. 
This is at variance with Ref. 1S ) , in which the inversion of the two levels was ob- 
tained by a particle-vibration coupling model with a non-standard Woods-Saxon 
parameterization for the mean field: a large diffuseness a=0.9 fm has been em- 
ployed and the spin-orbit parameter a\ s has been changed from 0.44 to 0.8 just 
for p-states. Thus, in such a modified Woods-Saxon potential, the lpi/2 and 2si/ 2 
states lie close enough so that the shift induced by phonons can reverse their order. 



We checked the above perturba- 
tive results for the | + state by di- 
agonalizing the Hamiltonian (2-1) in 
the five-dimensional model space with 
2si /2 , ld 5/2 02+, 2d 5/2 ®2+, ld 3/2 (g)2+, 
2d 3 / 2 <g>2 + configurations and found con- 
sistent results with the above calcula- 
tions. The wave function of the ground 
state is essentially given by 

1 + 

|- ) = 0.921|2s 1/2 ) + 0.356|ld 5/2 ®2+). 

(3-2) 

In the case of 24 core, the HF spec- 
trum with SIII interaction is shown in 
the first column of Fig. 3. We have 
then calculated + , 1~, 2 + , 3~ and 4 + 
vibrational states by the self-consistent 
RPA model. We first check the numer- 
ical accuracy of the model. The en- 
ergy and proton (neutron) moments M p 
(M n ) of the collective lowest 2 + state are 
obtained by three different methods for 
the RPA calculations. The energy and 
proton (neutron) moments are £'=3.90 
MeV and M p (M n )=1.8 (10) fm 2 , re- 
spectively, with unoccupied levels calcu- 
lated using the harmonic oscillator basis 

In the case of the coupling with the 
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Fig. 3. The neutron single-particle spectrum 
in 24 O core calculated with SIII-HF plus 
particle-vibration coupling. The energies 
are obtained, following the levels from left 
to right, by the HF potential, and includ- 
ing the effect of the 2^ state, of all 2 + 
states, and of all 2 + and 3~ states, respec- 
tively. The energies and B(EX) values of 
the excited states are calculated by the self- 
consistent HF+ RPA model. See the text 
for details. 



2f state, the ld 3 / 2 and 2si/ 2 states 
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have moderate shifts. Although the phonon is collective, there is substantial cancel- 
lation, in the case of the 2si/ 2 hole, between the diagrams (c) and (d) in Fig. 1 with 
the ld 5 / 2 and ld 3 / 2 intermediate state, respectively. In the case of the ld 3 / 2 particle 
the cancellation occurs between the two diagrams (a) and (b) with the ld 3/ / 2 and 
2si/2 intermediate states, respectively. There is not such cancellation for the ld 5 / 2 
state because the diagram (c) of Fig. 1, with intermediate 2si/ 2 state, gives a large 
effect and has negative sign, contrary to the general rule. The remaining quadrupole 
states have non- negligible effect although they are lying at higher energies. The 3~ 
states play a crucial role as far as the single-particle energy gap is concerned and 
deserves a special discussion. The large amount of low-energy (i.e., below 10 MeV) 
3~ strength explains the large absolute values of the energy shifts. In particular, two 
states at 6.07 MeV and at 7.33 MeV absorb 22% and 19% of the energy weighted 
sum rule strength, respectively. The fact that the 2sx/ 2 hole state is pushed down 
in energy is due to the absence of contributions from the diagram (c) of Fig. 1 be- 
cause of angular momentum selection rules, and also due to the strong contribution 
from the diagram (d) with the resonant-like lf?/ 2 intermediate state. These unusual 
effects are possible only in light nuclei where occupied states span only low values 
of the angular momentum, and were never found in systematic studies of heavier 
nuclei, for example, near 208 Pb. It is also interesting to observe experimentally the 
predicted collective 3~ states in 24 O to test the validity of the present theoretical 
study of the shell gaps. 

§4. Summary 

Effects of the collective modes on the shell structure have been investigated in 
the nuclei near the 10 Be and 24 core, at the neutron drip line. Energy shifts of the 
single-particle states are calculated by the particle-vibration coupling model taking 
into account the coupling to low-lying collective vibrational states. In the case of 
10 Be core, the coupling to the lowest 2 + state is found to be the most important effect 
to lower and raise the single-particle energies of 2si/ 2 and lpi/ 2 states, respectively, 
and leads to a much narrower gap between the two states than that of the Woods- 
Saxon potential. In the case of 24 core, the single-particle energies below the Fermi 
level behaves unusually in the calculations with the particle (hole)-vibration coupling 
compared with heavy nuclei near the 208 Pb core. The single-particle energy of the 
ld 5 / 2 state is lowered by the coupling to the collective 2 + states which gives an 
upward shift in heavy nuclei. Moreover the energy of 2s;l/ 2 state is also lowered by 
the coupling to 3~ states due to very specific effects found in the case of 24 O core, 
i.e., the blocking of the available configuration space for the particle(hole)-vibration 
couplings. 

Acknowledgments 

This work is supported in part by the Japanese Ministry of Education, Sci- 
ence, Sports and Culture by Grant-In- Aid for Scientific Research under the program 



6 



H. Sagawa , G. Cold and T. Suzuki 



number C(2) 12640284. 

References 

1) I. Tanihata et al., Phys. Lett. B287 (1992) 307, and references therein; 
I. Tanihata, J. Phys. G22 (1996) 157. 

2) P. G. Hansen, A. S. Jensen and B. Jonson, Ann. Rev. Nucl. Part. Sci. 45 (1995) 591. 

3) P. G. Hansen and B. Jonson, Europhys. Lett. 4 (1987) 409. 

4) H. Simon et al., Phys. Rev. Lett. 83 (1999) 496. 

5) H. Keller et al., Z. Phys. A348 (1994) 61. 

6) H. Iwasaki et al., Phys. Lett. B491 (2000) 8. H. Iwasaki et al., Phys. Lett. B481 (2000) 
7. 

7) D. Guillemaud-Moller et al., Nucl. Phys. A426 (1984) 37. 
T. Motobayashi et al., Phys. Lett. B346 (1995) 9. 

8) A. Ozawa et al., Phys. Rev. Lett. 84 (2000) 5493. 

9) X. Campi, H. Flocard, A. K. Kerman and S. Koonin, Nucl. Phys. A251 (1975) 193; 
I. Hamamoto, H. Sagawa and X. Z. Zhang, Phys. Rev. C53 (1996) 765; 

J. Dobaczewski et al., Phys. Rev. C53 (1996) 2809; 
J. Meng and P. Ring, Phys. Rev. Lett. 80 (1998) 460. 

10) A. Poves and J. Retamosa, Phys. Lett. B184 (1987) 311; Nucl. Phys. A571 (1994) 221; 
E. K. Warburton, A. Becker and and B. A. Brown, Phys. Rev. C41 (1990) 1147; 

S. E. Koonin, D. J. Dean and K. Langanke, Annu. Rev. Nucl. Part. Sci. 47 (1997) 463; 
Y. Utsuno, T. Otsuka, T. Mizusaki and M. Honma, Phys. Rev. C60 (1999) 054315. 

11) A. Bohr and B.R. Mottelson, Nuclear Structure, vol. II (Benjamin, New York, 1975). 

12) C. Mahaux, P.F. Bortignon, C.H. Dasso and R.A. Broglia, Phys. Rep. 120 (1985) 1, and 
references therein. 

13) G. Colo, Toshio Suzuki and H. Sagawa, Nucl. Phys. A695, 167 (2001). 

14) I.Hamamoto, H. Sagawa and X.Z.Zhang, Nucl. Phys. A626, 669 (1997). 

15) H. Sagawa and K. Asai, Phys. Rev. C63, 064310 (2001). 

16) When calculating the matrix elements of Hpv, the momentum dependent terms are 
dropped from v(r). 

17) V. Bernard and Nguyen Van Giai, Nucl. Phys. A348 (1980) 75. 

18) N. Vinh Mau, Nucl. Phys. A592 (1995) 33. 

19) F.M. Nunes, I.J. Thompson and R.C. Johnson, Nucl. Phys. A596 (1996) 171. 

20) A. Bohr and B.R. Mottelson, Nuclear Structure, vol. I (Benjamin, New York, 1969). 

21) F. Ajzenberg-Selove, Nucl. Phys. A490 (1988) 1. 

22) S. Raman et al., At. Data Nucl. Data Tables 36(1987) 1. 



